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Gauged B — L is a popular candidate for the origin of the conservation of R parity, i.e. R = 
(— ) 3 + i + 2 J j m supersymmetry, but it fails to forbid the effective dimension-five terms arising from 
the superfield combinations QQQL, u c u c d c e c , and u c d c d c N c , which allow the proton to decay. 
Changing it to B — XiL, where x e + x M + x T = 3 (with Xi 7^ 1) for the three families, would 
forbid these terms while still serving as a gauge origin of R parity. We show how this is achieved 
in two minimal models with realistic neutrino mass matrices, and discuss their phenomenological 
implications. 



Introduction 

In the Minimal Supersymmetric Standard Model 
(MSSM) of particle interactions, the imposition of R par- 
ity, i.e. R = (_) 3B + L + 2 i j where B, L, and j stand for 
baryon number, lepton number, and spin, respectively, 
serves two purposes. One is to establish a candidate for 
the dark-matter of the Universe, because the lightest su- 
persymmetric particle (LSP) is absolutely stable, being 
odd under R. The other is to forbid the otherwise allowed 
rcnormalizable superfield terms LH U , LLe c , LQd c , and 
u c d c d c , so that the proton does not decay as a result 
of these interactions. However, the higher-dimensional 
quadrilinear terms QQQL and u c u c d c e c are still allowed, 
giving rise thus to effective dimension-five terms in the 
Lagrangian which would also induce fast proton decay 
This is a serious problem for grand unification, even 
at a scale as high as 10 16 GeV. 

With the addition of three neutral singlet superfields 
N c , desirable for neutrino masses, gauged B — L becomes 
possible. In that case, R parity is better understood the- 
oretically as a discrete remnant of B — L breaking. How- 
ever, the offending quadrilinear terms remain, together 
with the new term u c d c d c N c . To get rid of these terms, 
we propose a simple solution. Instead of gauged B — L, we 
adopt a flavor-dependent variation, i.e. B — XiL, where 
the usual leptons still have L = 1, but x e ^ tT for the 
three families are not equal to 1. In fact, the idea that 
there is just one N c and that x e ^ tT — (3, 0, 0) or (0,3,0) 
or (0,0,3) has already been explored in the past [2|-l6|. 
If Xi 7^ 1, then all the unwanted quadrilinear terms are 
forbidden by gauged B — XiL so that proton stability is 
assured as shown for (0, 0, 3) case in Ref. Q. 

Here we assume that there are three N c , with each Xi 
nonzero. Unlike the previous (0, 0, 3) type of models, this 
opens up a new possibility that the spontaneous breaking 
of B — XiL by suitably chosen Higgs singlet superfields 
Si. 2 will result in an exact discrete residual symmetry 
which is just the usual R parity. In the following, we 
will discuss the conditions on Xj and construct two min- 



imal models with realistic neutrino mass matrices. We 
will examine their phenomenological constraints and the 
prognosis of their verification at the Large Hadron Col- 
lider (LHC). 



Model 

Our minimal model consists of the usual particle con- 
tent of the MSSM plus three singlet superfields Nf with 
L = — 1 and two singlet superfields Si^ with L = ±2. 
Under the gauged U(i)x symmetry of B — x^L, quarks 
have charges 1/3 and leptons have charges — Xi- The 
two usual Higgs doublet superfields have charges and 
the Higgs singlets Sijt have charges =F2^s- The vari- 
ous Xi and xs will be determined by the requirements of 
anomaly cancellation and a realistic neutrino mass ma- 
trix. 

The anomaly-free conditions for the addition of U(\)x 
to the standard SU(3)c x SU(2) L x U(l)y gauge group 
are easily written down. The [SU(3)c] 2 U(l)x and 
[U(l)x] 3 anomalies are automatically zero because of 
the vectorial nature of SU(3)c and U(l)x- The 
[gravity] 2 U(\)x anomaly also vanishes for the same rea- 
son. The [U(1)x] 2 U(1)y anomaly is zero because the 
sum of Y charges is zero separately for each family of 
quarks and for each family of leptons. The remaining two 



conditions, i.e. [SU(2)l] 
given respectively by 



'U(l)x and [U(l) Y } 2 U(l)x, are 



(3)(3)(l/3)-z e -z„-z T =0, (1) 

(3)(3)[2(l/6) 2 -(2/3) 2 -(-l/3) 2 ](l/3) 

+ [2(-l/2) 2 - (-l) 2 ](-x e - Xfl - x T ) = 0. (2) 

Both are satisfied if 



x e 



(3) 



The usual B — L is recovered for x e = x^ = x T = 1, and 
B — 3L T [3 is obtained for x e = i„ = and x T =3. 
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As Si t 2 acquire vacuum expectation values, U(l)x is 
broken. (x e , x^, x T ) as well as xs will be chosen in such 
a way that a residual symmetry of B — XiL will remain 
which is exactly R parity, and that realistic neutrino 
masses and mixing are obtained via the canonical see- 
saw mechanism. We note that since the Higgs doublets 
do not transform under U(l)x and Si ,2 do not transform 
under St/(3) c X SU(2) L x U(1) Y , the resulting Z and 
Z' x bosons do not mix. This avoids the stringent con- 
straint from precision electroweak measurements at the 
Z resonance. (See Ref. [7] and references therein.) 

To obtain R parity (equivalently, matter parity) as a 
total residual symmetry from the breaking of B — XiL 
using 5i 2i these new singlet superficlds (with L = ±2 
and B = 0) should satisfy \xs\ — 1/3. Since B = 1/3 
for quarks and L = 1 for leptons, we obtain the usual 
definition of R parity for all particles if 3x e .^. T are odd 
integers. For a detailed discussion about conditions to 
get a Zn out of a U(l) gauge symmetry, see Ref. [§] and 
references therein. 

Flavor-dependent U(l) models have been widely stud- 
ied to address many issues (e.g., see Ref. [§]). Our idea 
of having a particular discrete symmetry out of a flavor- 
dependent U(l) gauge symmetry can be viewed as a use- 
ful guide in constraining such models. 



Neutrino sector 

The general requirements discussed above would still 
allow an infinite number of possible models, until realistic 
neutrino masses and mixing are considered. 

If all X% QbTQ different for the three families of leptons, 
the charged lepton mass matrix and the Dirac mass ma- 
trix linking 2/, with N? are both constrained to be diag- 
onal. To obtain mixing among the three neutrinos, the 
3x3 Majorana matrix spanning N? must have enough 
nonzero entries. However, the only sources of such terms 
are S 1>2 N!fNg and ATjW| if Xj + x k = 0. With three dif- 
ferent i ejMjT and just ±x# to work with, this is clearly 
impossible. 

We now assume that x^ = x T , then some simple alge- 
bra will show that there are only two solutions (recalling 
that Si, 2 have L — ±2): 



(I) Xg — X e 

(II) x s = -Xf,, = (x e + Xfj)/2. 
Together with Eq. (j^J), this means that 



(4) 
(5) 



(I) x e!fi!r = (-3,3,3), x s = 3, (6) 
(II) Xe jM ,T = (9,-3,-3), x s = 3. (7) 



In model (I), the 3x3 Majorana mass matrix spanning 
Nf has all nonzero entries: N2N?, ^ are invariant mass 
terms, N^N^ comes from (S2), and T N^ T come from 



(Si). In model (II), the N£N£ entry is zero, but all others 
are nonzero: N° T N° T come from (S2), and N^N^ T from 
(Si). Both are general enough for obtaining a realistic 
neutrino mass matrix, with mixing among all three lepton 
families. 



Baryon triality 

As discussed above, the requirement of a realistic neu- 
trino mass matrix using only Si^ demands \xs\ — 3 in- 
stead of \xs\ — 1/3. This means that the total discrete 
symmetry from B — X{L is not just R parity any more. 
It has been extended to a larger symmetry Following 
the general arguments in Ref. [8j, we find that our total 
discrete symmetry is now Zq, which is a direct product 
of R parity and baryon triality (Zf, = R2 x B%). Un- 
der baryon triality (B3 = Z3) [10(, baryons transform as 
uj = exp(27ri/3), so that the proton is absolutely stable, 
being the lightest particle with that charge. This result 
came as a pleasant surprise, because it was not our inten- 
tion to construct a model which contains B%. It points 
to a possible deep connection between neutrino mass and 
proton stability. 
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anomaly- free U(l) gauge symmetry models i 
as well as discrete gauge symmetry models [1 
Especially, we note that when model (I) is shifted by 
some hypercharge, it can reach the equivalent form of a 
model in Ref. 11] . 



e — fi — r nonuniversality 

The salient prediction of our proposal is the existence 
of a new neutral gauge boson Z' x . It does not mix with 
the electroweak Z at tree level, but it couples to all 
quarks and leptons in a specified way. In particular, it 
breaks e — fi — t universality. Thus it may be important 
as a one-loop effect 0] in the precision measurements of 
Z — > £ + £~. However, this effect is proportional to a;f, 
and its contribution to nonuniversality is zero for model 
(I). As for model (II), the prediction is that T(Z — > e + e~) 
should be bigger than T(Z — > /i + /i~) and T(Z — > t + t _ ), 
and it is proportional to (81 
averages are [15 1 



9)g x . The present world 



T e = 83.91 ±0.12 MeV, 
r p = 83.99 ±0.18 MeV, 
T r = 84.08 ± 0.22 MeV. 



(8) 
(9) 
(10) 



After adding a kinematical correction of 0.19 MeV to F T , 
we find the deviation of T e from the average of T M and 
T r to be bounded at 95% C.L. by 



Ar e /r M , T < 0.002. 



(11) 
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Let r = M|, /M§, then the one-loop radiative correction 
to Z — > £ + £~ from Z' x exchange in model (II) is given 

by [Hi 



where 



I u,r °7T 



*a(r) 

2 

Lio 

6 Vl + r 



In 2 



1 + r 



(12) 



-- -2r- (2r + 3)lnr + 2(l + r) 2 x 



(13) 



10* 



10 J 



lO* 



10 1 
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In the above, Li 2 (x) = — (dt/t) ln(l — t) is the Spence 
function. For r > 1 (i.e. M§, > M§) which will be 
required by Tevatron data in any case (as shown below) , 
p 2 ~ r — 1 [ll/9 + (2/3) In?'] and the resulting numerical 
bound is 



9x/M 2 z , < 0.05 TeV" 



(14) 



to a good approximation. There is also an overall correc- 
tion to T(Z — > hadrons)/r(Z —> leptons), but that effect 
may be absorbed into the value of as in QCD. 



LEP2 contact interaction 



The agreement of the SM and the LEP2 data provides 
the indirect bounds on the Z' x mass. The strongest one 
for our models (I) and (II) , where /i and r have the same 
charges and e has the charge of opposite sign, comes from 



the e + e — > /i H 



channel with A vv = 16 TeV 



[17( | . (See Ref. [18| for some useful discussions.) 
The Z' x mass is bounded by 



M Z' X z ^Wx^\{hy V y 



(15) 



for sufficiently large Mz> compared to the LEP2 energy. 
It results 



(16) 
(17) 



M z > x £1.4 TeV in model (I), 
M Z ' X £ 2.3 TeV in model (II), 

for g x = 0.1. 

Phenomenology of Z' x 



The direct production of Z' x is possible at hadron col- 
liders through its coupling to quarks. Its decay to leptons 
is then a clear signature. At present, there is a Tevatron 
limit on the cross section a(pp — > Z' x — > e + e~) at a 
center-of-mass energy E cm — 1.96 TeV, based on an inte- 
grated luminosity of L — 2.5 fb -1 [lj|, and similarly for 
dimuons [20j . To compare against these results, we take 



FIG. 1: Tevatron bounds on Z' x mass in models (I) and (II). 



gx = 0.1 for definiteness. We assume Tz> — 0.01Mz> x 
for model (I) and ^z' x = 0.04:Mz' for model (II), which 
are approximately their true values if they only decay 
into particles of the Standard Model (SM), i.e. not their 
superpartners nor the additional singlets. 

For the numerical analysis, we use CompHEP/CalcHEP 
2lL l22j and the parton distribution functions of CTEQ6L 
Since the Tevatron bounds from dielectron and 



dimuon data are similar for the same coupling, we show 
only the dielectron result. We find Mz> £ 830 GeV for 
model (I) and Mz> x £ 940 GeV for model (II) as shown 
in Fig. [T] (Since IrcJ < \x e \ in model (II), the bound 
from the dimuon data is less stringent.) We note also 
that the e — /i — r nonuniversality constraint of Eq. (|14[) 
is easily satisfied. 

For the LHC discovery reach (with the design energy 
-E"cm = 14 TeV) through the dilepton Z' x resonance, we 
use cuts pt > 20 GeV, \rj\ < 2.4 (for each lepton) and 
\m im (l + (.-)\ < 3Tz x - SM background at the LHC with 
these cuts is negligible, and we just require 10 signal 
events to claim its discovery at the LHC for a fixed flavor. 

Figs.[2][a) and 2(b) show the LHC discovery reach us- 
ing dileptons for models (I) and (II) respectively. From 
the dilepton resonance only, model (I) cannot be dis- 
tinguished from the flavor-independent case of B — L. 
In model (I), Z' x will be revealed by both e + e~ and 
/i + channels at the same luminosity {L ~ 1 fb -1 
for M z > = 1-5 TeV). In model (II), the reso- 
nance will need a luminosity about an order of magni- 
tude larger than that for the e + e~ resonance because 
a(n+fi-)/a(e+e-) ~ (-3) 2 /(9) 2 . 

Since Higgs doublets have zero charges under B — XiL, 
the channels which require nonzero charges such as the 
6-lepton resonance discussed in Ref. [13] will be absent. 

In the presence of Z' x at the TeV scale, a (predomi- 
nantly right-handed) sneutrino as well as the usual neu- 
tralino can be a good LSP dark-matter candidate (25| . 
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FIG. 2: The LHC discovery reach for (a) model (I) and (b) model (II). The required luminosity is the same for both e + e~ 
and /i + A t_ resonances in model (I). The required luminosity for the e + e~ (solid line) resonance is about an order of magnitude 
smaller than that for the fJ. + ft~ (dashed line) resonance in model (II). 



One-loop charged lepton processes 

There are two sources of lepton flavor violation. One 
comes from explicit interactions linking jV c through £1,2; 
the other through the mismatch of lepton and slepton 
mass matrices (which is common to all supersymmetric 
models). Since we choose x e ^ x^ — x T , the latter ap- 
plies only to the fi — T sector, whereas the former applies 
to all leptons, but for one-loop processes such as fi — ¥ &y, 
they are negligible because of the smallness of neutrino 
masses. 

The muon anomalous magnetic moment (a^) does not 
violate lepton flavor, so it has a contribution from Z' x . 
i.e. El 



m. 



A ^ 



fx 



(18) 



for both models (I) and (II) . Using the Tevatron bounds 
on Mz 1 with gx = 0.1, we find this to be at most of 
order 10 -11 , which is negligible compared to the present 
experimental accuracy. On the other hand, since we are 
considering a supersymmetric model, there are one-loop 
contributions to et M from the neutralinos and charginos, 
which may account for the possible deviation in its mea- 
surement at the Brookhaven Alternating Gradient Syn- 
chrotron [271 ] . 



xs = 3 accomplishes three objectives, (i) The conven- 
tional R parity survives as an exact discrete symmetry, 
desirable for having the LSP as a good dark-matter can- 
didate, (ii) Realistic neutrino masses and mixing are ob- 
tained, (iii) An exact residual Z3 symmetry, i.e. baryon 
triality B3, emerges which makes the proton absolutely 
stable. 

The neutral gauge boson of this new U(l)x has large 
couplings to leptons over quarks: three or nine times 
larger than in the case of B — L. Using present Tevatron 
bounds with gx =0.1, we find Mz' to be greater than 
830 GeV in model (I), and 940 GeV in model (II). The 
indirect bounds from LEP2 are 1.4 TeV in model (I), and 
2.3 TeV in model (II). It should be accessible at the LHC, 
possibly at a very early stage. To have a background-free 
signal of 10 dilepton events at the LHC, the Z' x with 
Mz' x — 1.5 TeV in model (I) may be discovered for an 
integrated luminosity of only about 1 fb _1 . In model (II), 
because of the flavor-dependent charges, the event rate 
of the e + e~ resonance is predicted to be nine times that 
of n + /i~ . 

This work was supported by U. S. Department of En- 
ergy Grants No. DE-AC02-98CH10886 (HL) and No. 
DE-FG03-94ER40837 (EM). 



Conclusion 

We have proposed a new U{1) gauge symmetry B — 
XiL, with x e ^. T = (—3,3,3) [model (I)] or (9,-3,-3) 
[model (II)], in the context of supersymmetry with three 
neutral singlets Njj „ r . The spontaneous breaking of this 
U(l)x by the addition of singlets Si t 2 with L = ±2 and 



[1] S. Weinberg, Phys. Rev. D 26, 287 (1982). 
[2] E. Ma, Phys. Lett. B 433, 74 (1998). 
[3] E. Ma and U. Sarkar, Phys. Lett. B 439, 95 (1998). 
[4] E. Ma and D. P. Roy, Phys. Rev. D 58, 095005 (1998). 
[5] X. J. Bi, X. G. He, E. Ma and J. Zhang. laxXTv:0910.0771l 
[hep-ph]. 

[6] E. Salvioni, A. Strumia, G. Villadoro and F. Zwirner, 
lajXiv:091L1450l [hep-ph] . 



5 



[7] J. Erler, P. Langacker, S. Munir and E. R. Pena, JHEP 

0908, 017 (2009). 
[8] T. Hur, H. S. Lee and C. Luhn, JHEP 0901, 081 (2009). 
[9] P. Langacker and M. Plumacher, Phys. Rev. D 62, 

013006 (2000). 

[10] L. E. Ibanez and G. G. Ross, Nucl. Phys. B 368, 3 (1992). 
[11] K. S. Babu, I. Gogoladze and K. Wang, Phys. Lett. B 

570, 32 (2003). 
[12] H. S. Lee, Phys. Lett. B 663, 255 (2008). 
[13] K. S. Babu, I. Gogoladze and K. Wang, Nucl. Phys. B 

660, 322 (2003). 
[14] H. K. Dreiner, C. Luhn and M. Thormeier, Phys. Rev. D 

73, 075007 (2006). 
[15] [Particle Data Group] C. Amsler et al, Phys. Lett. B 

667, 1 (2008). 

[16] C. D. Carone and H. Murayama, Phys. Rev. Lett. 74, 
3122 (1995); Phys. Rev. D 52, 484 (1995). 

[17] [The LEP Collaborations ALEPH, DELPHI, L3, 
OPAL, the LEP Electroweak Working Group, 
the SLD Electroweak and Heavy Flavour Groups]. 
|arXiv:hep-ex/0312023| 



[18] M. S. Carena, A. Daleo, B. A. Dobrescu and 
T. M. P. Tait, Phys. Rev. D 70, 093009 (2004). 

[19] [CDF Collaboration] T. Aaltonen et al, Phys. Rev. Lett. 
102, 031801 (2009). 

[20] [CDF Collaboration] T. Aaltonen et al, Phys. Rev. Lett. 

102, 091805 (2009); 

[21] A. Pukhov et al.,|arXiv:h ep-ph/9908288 
[22] A. Pukhov, arXiv:hep -ph/0412191 

[23] J. Pumplin, U! RT Stump, J\ Huston, H. L. Lai, 
P. M. Nadolsky and W. K. Tung, JHEP 0207, 012 
(2002). 

[24] V. Barger, P. Langacker and H. S. Lee, Phys. Rev. Lett. 

103, 251802 (2009). 

[25] H. S. Lee, K. T. Matchev and S. Nasri, Phys. Rev. D 76, 
041302 (2007). 

[26] F. Jegerlehner and A. NyfTeler, Phys. Rept. 477, 1 
(2009). 

[27] [Muon g-2 Collaboration] G. W. Bennett et al, Phys. 
Rev. Lett. 92, 161802 (2004); Phys. Rev. D 73, 072003 
(2006). 



